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ABSTRACT 

This work extends our previous studies of two-photon annihilation into 
baryon-antibaryon pairs from spin-1/2 octet to spin-3/2 decuplet baryons. 
Our approach is based on perturbative QCD and treats baryons as quark- 
diquark systems. Using the same model parameters as in our previous 
work, supplemented by QCD sum-rule results for decuplet baryon wave 
functions, we are able to give absolute predictions for decuplet baryon 
cross sections without introducing new parameters. We find that the 
A ++ cross section is of the same order of magnitude as the proton cross 
section, well within experimental bounds. 
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1 Introduction 



The study of exclusive processes in quantum chromodynamics (QCD), where intact 
hadrons are explicitly measured in the final state, provides important insights into the 
mechanisms of confinement and into the dynamics of hadronic bound states [HI2|- 
Among the multitude of exclusive processes, two-photon annihilation into baryon- 
antibaryon pairs is particularly interesting, because it is one of the simplest calculable 
large-angle hadronic scattering reactions involving two hadrons. Therefore, 77 — > BB 
has recently received considerable experimental jH| and theoretical [U El IHj attention. 

In a recent paper [3] we have studied baryon pair production in two-photon col- 
lisions for baryons belonging to the lowest-lying flavor octet. In the present note we 
extend our work to reactions involving spin-3/2 decuplet baryons. Previously, two- 
photon annihilation into decuplet baryons has been studied in Refs. El EH E0 HD] 
within different frameworks with differing conclusions. Thus an experimental analysis 
could shed light on the relative importance of the underlying mechanisms considered 
here and in the aforementioned references. 

Our model is a modification of the perturbative hard-scattering picture (HSP) for 
exclusive processes [HI [TJ] . While the HSP is exactly valid only at asymptotically 
large momentum transfer, the interplay of perturbatively calculable with nonpertur- 
bative effects renders theoretical analyses quite intricate at energies where data are 
currently available. In order to parameterize such possible non-perturbative effects 
within a perturbative framework, an effective formalism was developed in Ref. |13j . 
where baryons are treated as quark-diquark systems. In the sequel this model has 
been successfully applied to a variety of exclusive reactions El EH EH1 HHj ITTj. 

In the following, we start with a brief review of the quark-diquark model. Then we 
go on to describe the new ingredients necessary for the study of processes involving 
decuplet baryons. In Sec. El we present and discuss model predictions with emphasis 
on the A cross sections, for which experimental upper bounds are available |18j . 
Following concluding remarks, supplementary analytical expressions for the scattering 
amplitudes are tabulated in the Appendix. 

2 Exclusive Reactions in the 
Quark-Diquark Picture 

Here we briefly summarize the modified hard-scattering formalism with diquarks, 
and elaborate on the aspects specific to the treatment of decuplet baryons. For a full 
account of all details we refer to our recent work [SJ [T7| . 
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2.1 Review of the Model 



As in the conventional hard-scattering picture [TTJ[T2j , an exclusive reaction amplitude 
M. is convolutively factorized into a process-dependent, perturbative hard-scattering 
amplitude T and process- independent, non-perturbative distribution amplitudes 
The latter are probability amplitudes for finding the pertinent valence Fock states, 
here quarks and diquarks, in the scattering hadrons. The amplitude for two-photon 
annihilation into a baryon-antibaryon pair is given by 

1 1 

M {x} {s,t) = Jdx 1 Jdy 1 ^i(x 1 )^t(y 1 )f {x} (x l ,y 1 ;s,t 



where Lorentz and color indices are suppressed for convenience. Furthermore, the 
dependence on renormalization and factorization scales is neglected since we are only 
interested in a rather restricted range of momentum transfer. The subscript {A} 
denotes all possible configurations of photon and baryon helicities. In the following 
we use the label B to denote spin- 1/2 octet baryons and B 10 to label spin-3/2 decuplet 
baryons. 

For the process 77 — > B 10 Bio, there are 19 independent helicity amplitudes, 
■M.\ B ;Ai,a 2 > where the Ab 10 , Ag io are the helicities of the outgoing baryon and 

antibaryon, respectively, and Ai, A2 label the helicities of the two photons. Only 13 
out of these 19 helicity amplitudes involve a zero or single flip of the hadronic helicity. 
Double flip amplitudes vanish in our approach. We use the following convention for 
the nonvanishing amplitudes: 



M_l Li 7 = M_i 3.! 

2'2> ' _ 2>2> > 

Al_1_1.il> 08 = -Mi -a-i i> 

2' 2' > __ 2' 2' ' 

^|-i;l,l> 09 = -Mi,-|;l,-l> 

•Ml.ijl,-!, 010 = AOf.1,1, (2) 



=__^,-_i 1,-1' 011 = -M_|, 

^6 = Mi, 1-1,1, 12 = ; M3 

013 = A^f,-!;!,-!- 



Other helicity configurations are related to these via parity and/or time reversal 
invariance. Our normalization of the amplitudes is such that the differential cross 
section for two-photon annihilation into decuplet baryons is given by 
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where the sum is over all possible helicity configurations {A}. 

In T consists of all possible tree diagrams that contribute to the elementary 
scattering process 77 — > qDqD. The momenta carried by quarks q and diquarks D are 
assumed to be collinear to those of their parent hadrons, B. The quark and antiquark 
carry momentum fractions x\ and y\ in the baryon and antibaryon, respectively, while 
the diquark and antidiquark carry momentum fractions x^ = 1 — x\ and y^ = 1 — yi, 
respectively. Since we assume that every baryonic constituent has a four-momentum 
xpb proportional to the four- momentum of its parent hadron p B [19], it acquires 
an effective mass xm-Q, where tub denotes the baryon mass. These effective masses 
are taken into account for all internal and external legs of the Feynman diagrams 
contributing to the hard-scattering amplitude T. The hard-scattering amplitude is 
then expanded in powers of the small parameter (itib/v^) u P to next-to-leading 
order, at fixed center-of-mass scattering angle 9. The result is reexpressed in terms of 
massless Mandelstam variables, s, t, and u which are obtained from the usual massive 
Mandelstam variables, s, t, u, again by expansion in (me/Vs). In the hard scattering 
diagrams, the composite nature of the diquarks is taken into account by diquark form 
factors. These are parameterized such that asymptotically the scaling behavior of the 
pure quark HSP emerges. 

The complete parameterization of the model, including form factors and octet- 
baryon wave functions can be found in [Sj. These parameters were fixed in ^3] by 
fitting elastic electron-nucleon scattering data. With the same set of parameters a 
variety of other processes has been computed, and the results have successfully met 
experimental comparison IHH E] • 



2.2 Decuplet Baryons 

The diquark model comprises spin-0 (scalar) and spin-1 (vector) diquarks. While 
both scalar (S) and vector (V) diquarks contribute to processes involving spin-1 /2 
octet baryons, the valence Fock states of spin-3/2 decuplet baryons consist only of 
quarks and vector diquarks. 

We recall that the valence Fock state of an octet baryon B with mass tub, mo- 
mentum pb, and helicity A can be described by the following quark-diquark wave 
function 

*b(pb,*, A) = /f$|(z)x|u(p B , A) + /B $ B (;r)x B 1 + p£ \ 75M(pB; A) (4) 

V3 V m B/ 

when transverse momenta of the constituents are neglected, x is the longitudinal 
momentum fraction of the quark, whereas the diquark carries the longitudinal mo- 
mentum fraction 1 — x. Analogously, the wave function of a decuplet baryon may be 
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written as 

KJpb 10 ,x,X) = /? 10 < (x)x? 10 m /1 (pb 10 ,A), (5) 
with the Rarita-Schwinger spinors j20j 

Vb 1 (p, A = ±3/2) 
u» (p,A = ±1/2) 

Recall that all Lorentz indices have been suppressed in Eq. ((TJ), the open index \i 
of the vector diquark polarization vector in (jlj) and © is contracted appropriately 
in the convolution integral (JTJ. Xd> Xd 10 (-^ = S,V) denote pertinent SU(3) quark- 
diquark flavor wave functions and $ B ) , $^ 10 represent the nonperturbative probability 
amplitudes for finding these constituents with momentum fractions x and 1 — x, 
respectively, in the (decuplet) baryon. These probability amplitudes are normalized 
such that 

i 

J dx$*(x) = 1, (7) 
o 

and analogously for $^ 10 . The constants /J, /^ 10 result from integrating out in- 
trinsic transverse momenta in the full wave function to produce Eqs. (0} and (jHJ), 
respectively. The numerical values of and f^ 10 are furthermore determined by the 
overall probability of finding the \qD)-sta.te in the baryon B or decuplet baryon Bi , 
respectively. 

For unbroken SU(6) spin-flavor symmetry octet- and decuplet baryon wave func- 
tions are related, specifically, $f = $^ = $Bio and y| = jb = /^icy^ Iri t h e 
actual parameterization of the diquark model the asymptotic SU(6) symmetry 
is systematically broken down to SU(3) flavor symmetry. Thus the above SU(6) re- 
lations are by no means satisfied, and $5 and $y as well as fg and fy have quite 
different values. Since SU(6) symmetry is thus already broken within the baryon 
octet we cannot use SU(6) symmetry for deriving quark-diquark wave functions of 
decuplet baryons. Instead, we will apply another strategy to fix $y 10 and fy 10 . 

The lowest moments of three-quark wave functions of octet and decuplet baryons 
are restricted by QCD sum rules [HI ISj- Model wave functions that satisfy the 
QCD sum-rule constraints (for a typical factorization scale of about 1 GeV) are 
very asymmetric in the longitudinal momentum fractions x,, % = 1,2,3 of the quarks 
for octet baryons and nearly symmetric (~ X1X2X3) for the As and the Q j^]. By 
regrouping terms in the three-quark wave function such that, for example, quarks 2 
and 3 are in a specific spin-flavor state and by integrating over one of the momentum 
fractions of the two quarks that build up this "diquark" we can convert the three- 
quark wave function into a quark-diquark wave function that nearly has the form 



£ M (±1) u(p,X = ±1/2) 



2A ( 



7s 



u(p,X). 



(6) 
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flU) or (JI3) for octet or decuplet baryons, respectively. For more information on this 
conversion we refer to jSj. The probability amplitudes $y and $y 10 for general three- 
quark wave functions are different in the cases of helicity-0 and helicity-1 V diquarks. 
For the octet and decuplet model wave functions that we employ this difference turns 
out to negligible. We then arrive at Eq. (jlj) or (JHJ), respectively. 

We apply the above procedure to the three-quark wave function of the A that 
has been proposed in Ref. j^j based on QCD sum-rule constraints. We obtain the 
following quark-diquark wave function for a A with helicity ±1/2 



1/2, 



X 



Nx(l 



X) 



[1 - 2.95x + 3.86x 2 ) exp <^ -b 2 



m„ 



x 



+ 



mi 



1 — x 



Analogous to the standard parameterization of the diquark model for octet baryons 
we have introduced an additional exponential factor that damps the end-point regions 
x — > 0, 1. Such an exponential factor results if the transverse momentum dependence 
of the full wave function, which is integrated over, is assumed to be of Gaussian form. 
The parameters b 2 = 0.248 GeV 2 , m q = 0.33 GeV, and my = 0.58 GeV are taken 
to be the same as for octet baryons. The normalization factor N is determined by 
Eq. (JJJ). The expression for $y'' A ' _3//2 (x) differs, in general, from $y'' A ' _1 ^ 2 (a;). How- 
ever, we refrain from quoting it here, because our explicit calculations show that the 
production of helicity-3/2 As is suppressed within the diquark model. 



The only remaining open parameter is now the normalization /- 



A,|A|=l/2 
V 



of the 



helicity-1 /2 A wave function. Since the normalization fy of the octet-baryon wave 
function was taken as a free parameter in the diquark model we normalize the A 
wave function relative to the proton wave function. This means that we convert 
the three-quark wave functions for proton and A into quark-diquark wave functions 
of the form (jlj) and (j3J), respectively, and consider the resulting ratio fy ^ 2 /f v . 
For the QCD sum-rule based wave functions of Refs. jU] and [21] this ratio becomes 
fy lM=1/2 /f v = 0.898. With f v = 127.7 MeV, the value obtained in a fit of elastic 
electron-nucleon scattering data [T3j, we thus find 



A,|A]=l/2 
V 



125.1 MeV. 



(9) 



This completes the parameterization of our model for decuplet baryons. For sake 
of completeness we quote the flavor wave functions entering (jSJ) for the differently 
charged As: 

U V{UU} 5 

V2dV {ud} + uV{ dd } \ I \Fi , 



A++ 
X V 

A+ 

Xv 



A 

Xv 



Xv 



dVs 



{dd} 



(10) 
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3 Results 



We list analytical results for the hard-scattering amplitudes T{a} contributing to 
77 — > B10B10 in the Appendix. These results have been checked via crossing relations 
[221 against the separately computed amplitudes for the crossed process, Compton 
scattering 7B 10 — > 7B 10 . Comparing the spinor structure of the decuplet baryon wave 
function (j3J) with the one for octet baryons (J3J), we find that the leading, non-flip, 
hard amplitudes for decuplet baryons with helicity 1/2 are related by a factor of 2 
to those for octet baryons. From the analytical expressions we also observe, that the 
hard-scattering amplitudes for octet baryons with helicity ±3/2 are suppressed by 
(9(m| io /s) or higher, even if these amplitudes conserve the hadronic helicity or flip it 
by one unit. The only 4-point contribution that is not suppressed enters the helicity 
amplitude 02- In the numerical calculations this contribution, however, turns out to 
be nearly negligible. 5-point functions with both photons attaching to the diquark 
do not contribute at all, since these are also suppressed by (9(m| io /s) or even higher. 

As a consequence of this observation, cross section ratios of different decuplet 
baryon channels can easily be estimated, provided that the corresponding probability 
amplitudes $y 10 are not too different. The cross-section ratios are then essentially 
determined by the corresponding charge-flavor factors C^' (see Eq. (fT3*jO and the 
wave function normalizations fy 10 . For the A-quartet $0 and fy are the same for all 
members due to isospin symmetry. From the flavor wave functions ([TT)|) the charge- 
flavor factors C c ( f 3) are seen to be 4/9, 3/9, 2/9, and 1/9 for the A++, A+, A , and 
A~, respectively. The cross section ratios become (approximately) 

a(A ++ ) : a(A+) : <r(A°) : a(A~) = 16 : 9 : 4 : 1 . (11) 

This is the first interesting prediction of the diquark model. In Fig. ^ we show 
the integrated cross sections (| cos(9cm)\ < 0.6, where 6cm is the center-of-mass 
scattering angle) for the A channels. The plot exhibits numerical predictions obtained 
with the standard parameterization of the diquark model [3] and the A wave function 
derived in Sec. 12.21 It confirms Eq. (fTTj) within 1 percent. 

This prediction is to be contrasted with the ratios 16 : 1 : : 1 that result if 
the photons couple to the total charge of the As. Also within the pure quark HSP 
the ratios for the A + and the A channels differ from ours. Within the pure quark 
HSP the cross section ratios for the different A channels are predicted to be cr(A ++ ) : 
cr(A+) : cr(A°) : a(A _ ) » 16 : 2 : 1/3 : 1 9J. Note that all the above predictions 
agree with our result for the cross section ratio a(A ++ ) : <j(A~) « 16 : 1. This result 
is also found in a more general QCD analysis However, yet another possible 

production mechanism via multi-pion intermediate states predicts <r(A ++ ) = cr(A~) 
and cr(A + ) = er(A°) [6j. An experimental determination of such cross section ratios 
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2 3 4 5 

W[GeV] 

Figure 1: Integrated cross sections for 77 — > A ++ A (solid line), A + A~ (dotted), 
A A (dashed), A~A + (dash-dotted line) (| cos(8cm)\ < 0.6) versus center-of-mass 
energy W = y/s predicted with the standard parameterization of the diquark model |H] 
and the A DA defined in the text (see Eqs. ©, ©, and ©). 



could therefore provide important clues on the underlying production mechanisms, 
especially because in ratios of cross sections for different A channels the sensitivity 
to the specific form of the A wave function should be greatly reduced. 

If we assume SU(3)-flavor symmetry, that is, if we take the same & v 10 and fy 10 
for all decuplet baryons, we are also able to give estimates for the pair production of 
strange decuplet baryons. Aside from appropriate phase space factors, SU(3) sym- 
metry implies 

a(A+) = <jpT+), 

a(A°) = a(£*°) = , 

o-(A-) = <r(E*-) = a(E*-) = a(n*~) . (12) 

However, since it is experimentally very difficult to measure pair-production cross sec- 
tions for decuplet baryons, we refrain from giving quantitative results for the strange 
decuplet baryons. We rather concentrate in the following on the A ++ channel which 
might have the best chance to be measured due to its comparably large cross section. 
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Figure 2: Integrated cross section for 77 — > A ++ A (|cos(0cAf)| < 0.6) versus 
W = y/s for the same A DA as in Fig El The solid line corresponds to the full 
diquark-model calculation. The contribution to the cross section that comes from the 
hadronic-helicity conserving amplitudes 4>i and 05 is represented by the dashed line. 
For comparison we also show the integrated cross section for 77 — ► pp (dotted line) 
calculated within the same model The shaded boxes indicate experimental upper 
bounds as obtained by the ARGUS collaboration [TH] . 



In Fig. El we show for comparison with the 77 — > A ++ A cross section the 
77 — y PP cross section that we have obtained with the same parameterization 
Surprisingly, we find that the A ++ cross section is of the same order of magnitude as 
the proton cross section. This prediction seems to be very stable against (reasonable) 
changes of the A wave function. With a A wave function that satisfies the SU(6) 
relations $y = $y and fy = v2/y we obtain, for example, a result which is only 
about 20% to 30% smaller *. Under the naive assumption that the photons couple 
directly to the charges of the baryons one would expect the A ++ cross section to be 
about 16 times larger than the proton cross section. 

*In a previous attempt to estimate a(A ++ )/a(p) within a diquark model a ratio of ~ 0.1 was 
found |BJ- This, however, was obtained with an incomplete version of the diquark model, where V 
diquarks were not taken fully into account and mass effects have been neglected. 
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From the viewpoint of the pure quark hard-scattering picture, the ratio of 77 — > 
A ++ A to the 77 — > pp cross section depends strongly on the choice of the proton 
wave function Not surprisingly, a result for the ratio comparable to ours is 
obtained with the QCD sum-rule wave functions of Refs. jH] and [21] for A and proton, 
respectively, which we have used in Sec. l2.2l to derive and normalize our quark-diquark 
wave function of the A. However, if the asymptotic wave function ~ 37x2X3 is taken 
for both the proton and the A, the cross section ratio a(A ++ )/cr(p) can be as large 
as 50 within the pure quark HSP [7|. On the other hand, soliton models involving 
multi-pion channels predict a much smaller ratio [21 E! , comparable to our findings. 
An experimental determination of the ratio a(A ++ )/cr(p) could therefore help to 
explore the importance of the various mechanisms that result in these quite different 
predictions. 

Unfortunately, it is very difficult experimentally to isolate the signal of the broad 
A ++ resonance from the background and to disentangle the A ++ and the A contri- 
butions in the 77 — > pp7r + 7r~ cross sections which are actually measured. Therefore 
only upper limits for the 77 — > A ++ A cross section have been extracted up to 
now by the ARGUS collaboration ^H]. As can be seen in Fig. [21 our results lie well 
below these upper limits. More recent attempts to constrain the 77 — > A ++ A cross 
section using the data taken by the L3 group are afflicted with the same problems, 
but a preliminary assessment indicates compatibility with the ARGUS results and 
our predictions [2SJ- A better chance to determine the cross section for A ++ A pair 
production would perhaps exist for the BABAR or BELLE experiments which enjoy 
a much higher luminosity. 

Finally, let us comment on the treatment of mass effects within our approach. 
Fig. [21 displays the effect of taking into account the finite A mass. As explained 
in Sec. 12.11 the A mass is taken into account in the hard-scattering amplitudes via 
an expansion in the small parameter (me/v^) where only the leading and next-to 
leading order terms are kept. As expected, mass correction terms do not contribute 
to the hadronic helicity-conserving amplitudes and 5 . Only the amplitudes that 
involve a single flip of the hadronic helicity, which vanish if masses are neglected, 
become nonzero due to the mass correction terms. The comparison of the solid and 
the dashed lines in Fig. [21 shows that these mass effects can be sizable in the few-GeV 
region. At W = 2.5 GeV the leading-order contributions provide only about 30% 
of the full cross section. This ratio increases, of course, with increasing energy and 
becomes roughly 70% at W = 5 GeV. 
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4 Concluding Remarks 



In this work we have computed 77 — > Bio Bio cross sections at intermediate momen- 
tum transfer for the case of spin-3/2 decuplet baryons Bi . We have employed a 
modification of the hard-scattering picture for exclusive reactions, where baryons are 
treated as quark-diquark systems, thereby effectively parameterizing nonperturbative 
contributions which are undoubtedly present at currently experimentally accessible 
energies. Using the same model parameters as in previous studies of other photon- 
induced reactions, and constraining the quark-diquark wave function of the A with the 
help of QCD sum-rule results, we are able to give absolute predictions for 77 — > AA 
without introducing new parameters. 

We find that the cross section for 77 — > A ++ A is of the same order of magnitude 
as the cross section for proton pair production, 77 — > pp. Furthermore, we observe 
that the pair production of decuplet baryons is almost completely determined within 
our model by those graphs where both photons couple to the quark line. This enables 
us to estimate production ratios for different decuplet-baryon channels independent 
of the choice of the wave function, provided that the wave functions are similar for all 
baryons within the decuplet. This is certainly the case for the A-quartet for which 
we predict the ratios a(A++) : <r(A+) : <r(A°) : a(A~) = 16 : 9 : 4 : 1. 

There are various other estimates of these cross section ratios in the literature, 
based on different viewpoints and production mechanisms, which differ in their pre- 
dictions from ours. It would therefore be necessary to compare to experimental anal- 
yses, in order to determine the relative importance of the considered production 
mechanisms, and to learn more about the degree of symmetry among constituents 
in decuplet-baryon distribution amplitudes. Such an experimental analysis should 
be quite feasible at a high-luminosity e + e~ collider. We therefore hope that our 
experimental colleagues will study this interesting problem in the near future. 

A Elementary Helicity Amplitudes 

for 77 — > qVqV 

There are 30 Feynman graphs that contribute to the hard-scattering amplitudes T 
for 77 — > qVqV. Their general structure is 

f {x} (i,u) = C®T?> V \t,u)FP + C$Tt D \i,uW + C7<? Tf V) (U)if \ (13) 

where are the appropriate charge-flavor factors. The subscript i = 1, ... ,13 
labels the helicity-combinations according to Eq. (0) . Each n-point contribution T <n,V '' ) 
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is found from a separately gauge-invariant set of Feynman diagrams, where (n — 2) 

( n y) 

gauge bosons couple to the diquark. The T ' are multiplied with the appropriate 
diquark form factors Fy \ parameterizing the composite nature of diquarks. For 
further details we refer to jSJ. 

The analytical results for T- n ' V are presented in the following. For their calcula- 



tion we employed the algebraic computer program Mathematica [24J with the package 
FeynCalc [22] • We do not list n-point contributions that are suppressed by at least 
O(m\ 10 / s), since they are neglected in our numerical calculations. Those include, 

for example, all 5-point functions Tf' V ^ and all amplitudes with (anti)baryon helic- 
ity ±|. Note that the parameterization of the form factors Fy 1 ^ for 4- and 5-point 

(3) 

functions provides additional inverse powers of s as compared to F v . We abbreviate 
C = (47r) 2 Cira a s , where Cp — § is the color factor, and a denotes the fine structure 
constant a ~ 1/137. Ky is the anomalous magnetic moment of the vector diquark. 



T x [t, u) 



T 2 (t, u) 
T 4 (t, u) 



4c- 



u 



3 m^ w Vut \XiVi + x 2 y2, 

2^ 1 V§sx x +y x 
3 m Bl0 ui x x y x 
c 2 1 If 



3m Bl0 \rsutx\X2y\y2 [ 
-Ky (2x x — S)y2t 2 + (2y x — 3)x2U 2 — Ax x y x ui 



+ 



+ 



,%2 



+ 1/2) (i 2 x x y 2 + u 2 x 2 y x - 2x x y 2 ut) - y 2 V 



x 2 u 



Tp(3,V) a r v 

T 5 (t, u) 
T 6 (i, u) 



(i,u) 



T x (u,t) 



-C 



2 1 



s s 



3 m Bl0 wt 
2 ^y(l ~ K v 



. ^x x yl + y x xl x x +y\ 

1 + « v ) h «y 



xix 2 y x y 2 



x 2 y 2 



mi 



xix\y x y\' 



The hard-scattering amplitudes for Compton scattering off decuplet baryons are 
related to the amplitudes listed above via crossing [22] • The corresponding elementary 
helicity amplitudes, 'yqD — > 7<?-D have been computed separately as a check. They 
can be obtained from the authors upon request. 
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